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Abstract. The balance function, which measures the correlation between opposite sign charge 
pairs, is sensitive to the mechanisms of charge formation and the subsequent relative diffusion 
of the balancing charges. The study of the balance function can provide information about 
charge creation time as well as the subsequent collective behavior of particles. In this paper, we 
present a reaction-plane-dependent balance function study for Au+Au collisions at y / %N = 200 
GeV and compare with results from recent three particle correlation measurements. We also 
report balance functions for relative pseudorapidity (A77), relative rapidity {Ay), and relative 
azimuthal angle (A<jj) from the recent RHIC beam energy scan data. 



1. Introduction 

The study of correlations and fluctuations can provide information about the properties of 
the hot and dense matter created in relativistic heavy-ion collisions, one such observable, the 
balance functions, which measures the correlation between the opposite sign charge pairs, are 
sensitive to the mechanisms of charge formation and the subsequent relative diffusion of the 
balancing charges [ij. Due to conservation laws like electric charge conservation, particles and 
their anti-particles are pair produced and correlated initially in coordinate space, if a delayed 
hadronization occurs, the lower temperature and less expansion and diffusion will result in a 
narrower charge balance function. 

The reaction-plane-dependent balance function in azimuthal angle can be written as 



r>(± a ja l r A + _(0,A0)-A ++ (0,A0) A_ + (0,A0)-A__(0,A0) 1 

Here N + ^((j)) is the total number of positive (negative) particles that have an azimuthal angle 

(f) with respect to the event plane, A^ (<p,A(p) represents the total number of pairs summed 

over all events where the first (positive) particle has an azimuthal angle 4> with respect to event 
plane and the second (negative) particle has a relative azimuthal angle Acj) with respect to the 
first particle. Similarly we can express A ++ ((fi, A(ft), A |_(c/>, A<p) and A ((/>, A<p). 

On the other hand, it has been discussed recently that the hot and dense matter created 
in heavy ion collision may form metastable domains where the parity is locally violated, this 
possible local parity violation [2] coupled with strong magnetic field produced by passing nuclei 
in such a collision could cause a charge separation across the reaction plane in non-central 



collisions called the chiral magnetic effect (CME) [3jlU[5]. ^ ne observable proposed to measure 
the Chiral Magnetic Effect is the three point correlator [6] , where <f> is the azimuthal angle of a 
particle and ipRp is the reaction plane angle. 



7a,/3 =< cos((j) a + 4>fS - 2tp RP ) > . (2) 

2. Experimental Method 

The data used in this analysis is from Au+Au collisions at -^/snn = 200 , 62.4, 39, 11.5, and 7.7 
GeV taken by the STAR experiment. The Time Projection Chamber (TPC) was used as the 
main detector for charged particle tracking. All tracks were required to have a distance of closest 
approach (DCA) to the measured event vertex of less than 3 cm. For all charged particles, a 
transverse momentum cut of 0.2 < pt < 2.0 GeV/c was applied while for identified particles, 
we used a pt cut of 0.2 < p t < 0.6 GeV/c. The pseudorapidity cut used is \rj\ < 1.0. Particle 
identification was done by selecting particles whose specific energy losses (dE/dx) were within 
two standard deviations of the energy-loss predictions for a given particle type and momentum. 
In addition, electrons were excluded from the analysis for all cases by specific energy loss inside 
the TPC. To correct for differences between the acceptances for positive and negative particles, 
a mixed event subtraction is used for all balance function results throughout this paper. 

To determine the event plane angle, the second Fourier harmonic in the azimuthal angle 
distribution is used. We used all tracks from the TPC while a pt-weight method was used to 
maximize the event plane resolution. Also the (j) weight method is applied to flatten the event 
plane distribution. 

3. Results 

3.1. Event- Plane- Dependent Balance Function at y / SNN = 200 GeV 

Figure [1] shows eft = 0° (in-plane) , (j> = 45° , and eft = 90° (out-of-plane) balance function for 
40-50% centrality only. The in-plane balance function is narrower than the out-of-plane balance 
function, which is caused by the stronger collective flow in-plane: charge pairs are created closely 
in space and time due to conservation laws, this correlation would remain to the final stage if 
they are not diffused, which gives a narrower balance function. Figure [T] also shows the <j> = 45° 
balance function, which is asymmetric and peaked at negative A(f>. Because of the strong elliptic 
flow established at non-central collision, if charge pairs are emitted with an angle respect to the 
reaction plane, the correlation would be stronger on the in-plane side compare to out-of-plane 
side. Also shown are the blast- wave model calculations [7], we can see that for all three cases 
here, the blast-wave model agree well with experimental data. 

To quantify the collective flow effect on balance function, we also study the weighted average 
cosine, Cb(cft), and sine, Sb(4>), extracted from the balance functions. 



= ~TJ\ i dA0£(0, A(/>)cos(A0), 
ZbW) J 

= 7TT / dA05(0,A0)sin(A0), 
ZbW) J 

z b (<j>) = [ dA0B(0,A0). (3) 



Cfc(</>) represents the width of balance function. If charges are created at the same point and 
did not diffuse due to strong collective flow, (%(</>) would be close to unity. s^((j)) is an odd 
function of A(j), so it can quantify the asymmetry of balance function. Figure [2] shows c&(0) 



and Sb{4>) for Au+Au collisions at y / SNN = 200 GeV. In the figure, Cb(4>) is closer to unity in 
the 0-5% centrality bin, which is due to a stronger collective flow in central collisions, while 
in mid-peripheral and peripheral collisions, Cb{4>) shows a difference between the in-plane and 
out-of-plane balance functions, which is caused by stronger elliptic flow in-plane. For Sb((f>), 
it reaches maximum at </>=135°, 315° and minimum at </>=45° , 225° , which demonstrates that 
charged pairs are more correlated on the in-plane side when emitted with an angle 4> with respect 
to the reaction plane. 
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Figure 1. (Color online) The balance 
function for = 0° (in-plane), eft = 45°, and 
<p = 90° (out-of-plane) particles pairs. The 
40-50% centrality bin is shown. The points 
are from the experimental data (not corrected 
for event plane resolution), while solid lines 
represent the blast-wave model calculations 



Figure 2. (Color online) The weighted 
average cosine and sine of balance function, 
four centralities are shown here. The points 
are from the experimental data (not corrected 
for event plane resolution), while solid lines 
are from the blast- wave model 171. 



Figure [2] also shows a comparison with the blast- wave model [7J. The blast- wave model 
includes a breakup temperature Tki n , the maximum collective velocities in the in-plane and 
out-of-plane directions, the spatial anisotropy of the elliptic shape by fitting STAR published v 2 
and spectra data [8], this model also assume local charge conservation and initial separation of 
balancing charges at freeze-out by fitting experimental results [9]. The difference between data 
and the blast-wave model could due to the finite event plane resolution for the data. 

The difference between the same-sign and opposite-sign three point correlator can be 
expressed as [7] 



1 2 

lv = 2^ 7+ " ~ 7++ ~ 1 — ) = M^ 2< Cb ^ > +V2c ~ V2s ^ ^ 

where 

V2c =< C b ((j)) cos(2cf)) > -v 2 < c b (4>) >, 
V2s =< s b ((f>) sin(20) >, 

and the bracket represents 
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Figure 3. (Color online) parity observable Figure 4. (Color online) The balance function 



7p scaled by experimental multiplicity 



in terms of Arj for all charged particle pairs from 
Au + Au collisions at snn = 39 GeV for nine 
centrality bins. 



In this equation, V2 {cb{4>)) will be positive if there are more charge pairs in-plane than out- 
of-plane, V2 C will be positive if the charge pairs are more correlated in-plane than out-of-plane, 
while V2s will be negative if the charge pairs are more correlated on the in-plane side. 

Figure [3] shows the parity observable calculate from balance function as well as its three 
components. All data points are corrected for the event-plane resolution here. To compare with 
previous results, we also plot the jp from STAR published data [TO] scaled by the measured 
uncorrected multiplicity in the same plot. Mathematically, the balance function result should 
equal the one from 7p and they do agree well. We can also see that, A thermal blast-wave model 
[7] incorporating local charge conservation and flow reproduces most of the signal. 



3.2. Balance Function at y / SNN = 39 GeV 

Previously, STAR has published measurements of the balance function from Au + Au, d + Au, 
and p + p collisions at y / SNN = 200 GeV [9j. In that paper, a simplified(Eq.([5]), no reaction- 
plane-dependence) balance function was used to measure the charge pairs correlation in terms 
of relative pseudorapidity, Arj, relative rapidity, Ay, relative azimuthal angle, Ac/>, and invariant 
relative momentum, qi nv 
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It has been shown that the balance function in terms of A77 and A(f> for all charged particles 
and in terms of Ay for charged pions narrow in central Au + Au collisions. Here we perform a 
similar measurement from Au+Au collisions at ^/snn = 39 GeV. 

Figure^] shows the balance function in terms of A77 for all charged particles for nine centrality 
bins. The balance function for central collisions is much more narrower compare to the one from 
peripheral collisions. The balance function for shuffled events is much wider than the balance 
functions from the measured data. The shuffled events are created by randomly shuffling the 
charges of the particles in each event. By doing this, all charge- momentum correlations are 
removed in shuffled events, which could produce the widest balance function within the same 
experimental acceptance. 

The balance function in terms of A</> would be more sensitive to transverse flow effect. FigureO 
shows the balance function in terms of A(p for all charged particles for nine centrality bins. The 



data shows a peak at around A<p = in central collision while it is almost flat in peripheral 
collisions. This is due to strong radial flow in central collision, when charge pairs are created, 
they are strongly correlated in space and time due to conservation laws, if a flow of perfect liquid 
established in the system, charges would remain correlated throughout the whole process which 
would provide this observed peak at A<j) = 0. Also due to STAR's uniform rj acceptance, the 
balance function for shuffled events are always flat for all centralities. 

Study of identified particle's balance function could reveal the charge production mechanisms 
of different particle species. Figure show the balance functions for identified charged-pion pairs 
for nine centrality bins in terms of the relative rapidity. The balance function for experimental 
data gets narrower in central collisions, while balance functions for shuffled events are much more 
wider. The dip of balance function for Ay < 0.2 should due to interpair correlations (HBT, etc) 
according to model predictions 
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Figure 5. The balance function in terms 
of A(fi for all charged particle pairs from Au 
+ Au collisions at ^/snn = 39 GeV for nine 
centrality bins. 



Figure 6. (Color online) The balance function 
in terms of Ay for identfied pion pairs from Au 
+ Au collisions at y / SNN = 39 GeV for nine 
centrality bins. 



3.3. Beam Energy Dependence of the Balance Function Widths 

To quantify the correlation measured from balance function, we calculate the weighted averages 
of < Ar/ >, < A(f> >, and < Ay > using Eq.Q. To remove the interpair correlation, the 
weighted average is calculated within 0.1 < Aij < 2.0 for < Ar] > and 0.2 < Ay < 2.0 for 
< Ay > 
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The left panel of Figure [7] shows the balance function widths for all charged particle in 
terms of the number of participating nucleons at ^/snn = 200 , 62.4, 39, 11.5, and 7.7 GeV. 
For both Ar] and A(f), the balance function widths narrow in central collisions for all energies 
presented here. The upper right panel of Figure [7] shows the energy dependence of the balance 
function widths of Ar] for most central events(0-5%). The observed < Ar] > from data shows 
a smooth decrease with increasing beam energy, which is consistent with the effect that radial 
flow is stronger at higher collision energies. A UrQMD [12J model calculation incorporating 
STAR acceptance filter show the similar trend but over predict the value, which is due to early 
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Figure 7. The balance function widths in terms of Arj, Acp for all charged particle pairs from 
Au + Au collisions. Left panel shows centrality dependence, while the right panel shows beam 
energy dependence of most central evens (0-5%). 



hadronization time and less flow in UrQMD. An early hadronization time would cause a wider 
balance function, while relative diffusion and interaction after freeze-out could further broaden 
it. The lower right panel of Figure [7] shows the same thing for < >. The experimental data 
shows a decreasing balance function widths with increasing beam energy, while UrQMD's width 
shows a little increase with increasing beam energy. This could due to final state interaction of 
produced particles in the UrQMD model. 

Figure [8] shows widths of the balance function for identified charged pions and identified 
charged kaons for both centrality dependence and beam energy dependence. For identified 
pions, the measured balance function widths get narrower in both central collisions and higher 
energies, the UrQMD model shows a similar energy dependence but over predicts the signal. In 
contrast, the widths of the measured balance function for identified kaons show little centrality or 
beam energy dependence, while it is consistent with the predictions of UrQMD model. This lack 
of centrality and beam energy dependence could indicate the different production mechanism 
of kaons: they are produced mainly at the beginning of the collision rather than during a later 
hadronization stage. However, the balance function for kaons is also affected by the decay 
<p K + + AT~, a further study of balance function for qi nv could help to distinguish the 
contribution from cp decay. 

4. Summary 

In this paper, we have presented new results from the reaction-plane-dependent balance function, 
the reaction-plane-dependent balance function analysis gives the same difference between the 
like-sign and unlike-sign charge dependent azimuthal correlations as the three point correlator 
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Figure 8. (Color online) The balance function widths in terms of Ay for identified pions and 
kaons from Au + Au collisions. Left panel shows centrality dependence, while the right panel 
shows beam energy dependence of most central evens (0-5%). 



results published by STAR. A thermal blast-wave model incorporating local charge conservation 
and flow can reproduce most of the difference between like- and unlike-sign charge-dependent 
azimuthal correlation. 

We also measured the none event-plane-dependent balance function for Au+Au collisions at 
y / «NN = 200 , 62.4, 39, 11.5, and 7.7 GeV for all charged particles, identified charged pions, 
and identified charged kaons. We find that the balance functions in terms of Aij and A<j> for 
all charged particles narrow both in central collisions and higher collision energies, which is 
consistent with the delayed hadronization picture. The UrQMD model reproduced the trend for 
energy dependence, but over predicted the value of balance function widths. For identified kaons, 
the balance function widths show little centrality or beam energy dependence, which indicate 
that kaons might be created early in the collision. However, a detail study using balance function 
for 

qinv is necessary to remove the effect from <fi decay. 
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